Scaling Limits for the 2D Metal-Insulator Transition at B = in Si-MOSFETs 
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We have reexamined data on the possible two dimensional metal-insulator transition at B = in 
Si-MOSFETs using a nonlinear regression method to extract all scaling parameters in a single fit. 
By keeping track of the magnitude of errors in the data we can use the normalized mean square 
deviation \ 2 °f the fit as a quantitative measure of how well the data is compatible with scaling. 
We have used this method to study electric field scaling in three different samples. We find rather 
good agreement of the data with scaling in individual fits, i.e. x 2 ' s °f about f , but also rather large 
variations in the fits depending on how cut-offs are introduced in the data. In particular, we report 
how fitted parameters vary when we cut away data that are either far from the critical point or 
at low excitation power, where temperature effects presumably dominate. In this way we find the 
critical i?-field exponent /3 to vary from about 3.2 to 3.9 with considerably smaller statistical error 
estimates in each fit. 



I. INTRODUCTION 

Several years ago the observation of a conductivity 
phase transition in the two-dimensional electron sys- 
tem (2DES) in Si metal-oxide-semiconductor field-effect- 
transistors (MOSFETs) was reportedtl. This observation 
was unexpected and controversial because it was contrary 
to thej-svidely accepted scaling theory for a noninteracjt 
2DES0. Further expexknental studies on MOSFETsf 
SiGe heterostructuresQu, and GaAs heterostructurescl~E£l 
indicate that a conductivity transition is possible as a 
function of carrier density in a 2DES at B = 0. The ob- 
served conductivity transitions at B = bear a marked 
similarity to transitions between states with different Hall 
resistances in the quantum Hall effect (QHE)Btjii3. This 
unexpected behavior and similarity to QHE transitions 
has also sparked considerable theoretical interestEZrEZI. 
Furthermore, these parallels make it advisable to con- 
sider developments in the QHE when studying the B = 
conductivity transition in 2DES. In particular, a recent 
publication by Shahar et a/Hreports that the behavior of 
the resistivity as a function of temperature and electron 
density, p(T,n s ), is inconsistent with scaling at tempera- 
tures above 0.5 K. Thus, they question the use of scaling 
analysis for the study of the QHE in general. We feel 
that a serious consideration of the appropriateness of the 
scaling analysis for the B = conductivity transition 
in 2DES is warranted. We have data available that can 
address this question using a set of data manipulation 
routines in Mathematica™ especially developed for this 
purpose. Below we briefly describe the data used in this 
study and the data manipulation routines that allow us 
to consistently propagate errors and eventually place con- 
fidence limits on the applicability of this scaling analysis 
to the case of the B = conductivity transition in 2DES. 
We feel that these methods have applicability beyond the 



particular system to the QHE and other quantum phase 
transitional. 

Our aim has been to develop an objective method that 
gives a quantitative answer to the question of how well 
the data actually scales. We achieve this by analyzing 
our ability to fit the data to a general scaling model that 
assumes as little as possible about the functional form 
of the scaling function. We use the variance-weighted 
mean-square deviation from this fit, x 2 : as a quantita- 
tive measure of the 'goodness' of the fit. Such a prop- 
erly weighted and normalized x 2 close to one indicates 
good agreement of the data with the model. An ob- 
jective estimate of y 2 requires a quantitative knowledge 
of experimental uncertainties. Therefore, great care has 
been taken to estimate the magnitude of these sources 
of uncertainty and to propagate them consistently to the 
final fit. 

Using this method to investigate E'-field scaling of the 
lowest temperature data available, we find rather good 
agreement with scaling in the best sample with a % 2 
varying from 0.4 to 1.9 depending on our selection of 
the data. We find that \ 2 gets smaller as more data at 
low power is discarded. Excluding more data at lower 
powers also tends to increase the consistency among the 
different samples. 



II. DESCRIPTION OF THE DATA INCLUDING 
SOURCES OF ERROR 

The data used in this study was amassed at the Uni- 
versity of Oklahoma while studying the temperature 
T, electron density n Sl and excitation dependence of 
the resistivity p of a number of high-mobility (/^ max > 
25'000 cm 2 /Vs) metal-oxide-semiconductor field-effect 
transistors (MOSFETs). These devices are described in 



detail elsewhere.l 

This data includes large sequences of currents / with 
their corresponding probe voltages V recorded as a func- 
tion of gate voltage V g oc n s and T using a four-terminal 
technique. The resulting /-Vcurves at low temperatures 
are highly nonlinear with distinctly different behaviors 
above and below a critical density n c or resistivity p c . 
In this article we consider the nonlinear /-Vcurves taken 
at the lowest temperature (~ 0.2 K) and present a scal- 
ing analysis of the resistivity as a function of the applied 
electric field, E. This E is taken to be the measured V 
divided by the distance between the respective probes. 

We find three sources of randomness that we assume 
are independent: 

• The uncertainty in the measured probe voltage 
characterized by a root mean square (RMS) un- 
certainty ay. 

• The uncertainty in the measured excitation cur- 
rent, 07, which we find is negligible. 

• The uncertainty in n s due to the applied V g simi- 
larly characterized by a n . 

For these samples at this T we also have a source of sjffi-i 
tematic error that we associate with the well-knownE2l 
difficulty in producing ohmic contacts to a 2DES. The 
non-ohmic contacts caused the measured I-Vcmves to 
be asymmetric with respect to I — and to show a sharp 
anomaly centered at / = 0. By carefully accounting for 
all the voltage drops in the measurement circuit, we were 
able to determine the contribution due to the contacts. 
We found that a percentage of this contact potential (pro- 
portional to the resistivity of the sample) added to the 
measured V essentially eliminated the observable effects 
of the contacts, see Fig. El Furthermore, we found that 
these contact potentials did not effect our analysis be- 
cause they occurred at very low excitation currents where 
the physics of the system is still dominated by the finite 
T. The results presented here are independent of these 
contact anomalies and our treatment of them. 

In order to measure the effective ay for our data, we 
produced a modified data set that was equally spaced 
in excitation I by using a spline interpolation scheme. 
The Fourier transform of the data was then examined. 
The transformed data is consistent with a low frequency 
signal superimposed on uniform white noise of constant 
amplitude. This white noise amplitude varies from about 
2 pV in the most metallic curves up to 90/iV in the most 
insulating ones. We used this amplitude as a measure of 
ay for each particular V g in the following analysis. 

The uncertainty in the electron density, a n , is more 
difficult to estimate. To calibrate n s to V g , we sweept 
the gate voltage at a constant magnetic field to ob- 
tain p X x(Vg). Minima in p xx that correspond to known 
integral filling factors, and therefore known n s , were 
then determined. The linear fit to the resulting (n s , V g ) 
pairs that gives V g (n s ) with a statistical error of a n = 
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FIG. 1. 
E. The upper curve taken at electron density n s 
cm -2 is in the insulating region, whereas the lower curve at 
n s — 1.05 * 10 11 cm -2 shows metallic behavior. The grey 
curves are corrected for the contact anomaly by adding a small 
fraction of the contact potential drop to the measured signal. 

1 * 10 9 cm~ 2 . This corresponds to a resolution of 0.8 
millivolt in V g . This determination is not obviously a 
measure of the density fluctuations for constant V g . Fur- 
thermore, an effective gradient in the gate voltage will 
be imposed on the 2DES by the source-drain excitation. 
We made every attempt to make this potential drop sym- 
metric with respect to the V probes, but the n s gradient 
could not be fully eliminated. Thus, we estimate the den- 
sity uncertainty, a ni by taking ay ~ lmV. This estimate 
is compatible with our observations of the precision and 
reproducibility of producing n s with V g in our samples 
as evidenced by a repeatability of a given p xx (V g ) where 
Pxx(n s ) is particularly steep. n 

Qualitatively, as has been reported previouslyp two 
different nonlinear V(I) behaviors can be seen. At low 
densities the resistivity has a maximum at zero electric 
field and decreases for stronger fields. This is similar to 
a normal insulator where a larger electric field can excite 
more carriers and hence increase the conductivity. At 
high electron densities the opposite behavior is seen; the 
resistivity has a minimum at E = and increases for 
stronger fields. Between these behaviors, a linear V(I) 
is observed with a constant p, see Fig. 0. This is similar 
to the nonlinear region of a superconductor. However, 
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FIG. 2. Nonlinear resistivities V/I for electron densities 
varying from 0.90 to 1.14 * 10 11 cm -2 . Here a low power 
cut-off P c — 50pW is introduced, discarding all data taken 
at a total power of 50pW or less dissipated in the sample. 
Other possible low power cuts are indicated in grey. Here 
no contact correction for the anomalies at the center of the 
curves is used. 



the nonlinearities seen here are anomalously strong, they 
occur at very small fields, and they are unexplained by a 
clear microscopic model. 

Because the I-V curves are nonlinear, there is a dif- 
ference between the normal resistivity given by the ratio 
V/I, and the differential resistivity, 4j. We have ana- 
lyzed both, and they have qualitatively similar behaviors. 
The differential resistivity has the advantage of being in- 
sensitive to experimental voltage and current offsets that 
must be considered when calculating ratios. The ^j is 
also more sensitive both to structure in the I-Vand to 
noise. We generate differential resistivities by doing a 
n-point differentiating convolution of the 1-Vda.ta,; 



Pi =^2 C j( V i+3 - Y i- 



(1) 



where the data points, {V^}, are made equidistant in cur- 
rent excitation by spline interpolation. Choosing differ- 
ent Super-Lancos differentiation filters, {c n }, allows us 
to optimize the data smoothing to reduce noise while re- 
taining the signal. We then propage ay in the original 
data through the convolution to obtain the variance in 
the resistance, o p . 



III. SCALING MODEL 

The question at issue is whether or not a quantum 
(T = 0) phase transition between two different conduc- 
tivity states takes place at B = in these systems. As 
mentioned in the introduction this is still highly debated 
in spite of large experimental and theoretical efforts. Our 
aim here is to develop a more rigorous and objective way 
to determine the consistency of experimental data with 
such a phase transition. 



A quantum phase transition occurs at zero temper- 
ature whereas all conductivity measurements are taken 
at both finite bath temperature and nonzero excitation 
voltage to drive a current through the sample. Theo- 
retical arguments suggest that an underlying quantum 
phase transition should give rise to. scaling at finite tem- 
peratures and excitation £ , -fields,E3 i.e. that the resistiv- 
ity of the sample should have a specific temperature and 
_E-field dependence; 



p(n s ,T,E) = p{x,y) = f 



£1/0 ' T l/uz 



(2) 



where 8 — " s " c is the relative deviation of the electron 
density from a critical density n c , and v, z and fi are crit- 
ical exponents associated with the phase transition. As- 
suming the simplest E'-field heating mechanism, j3 would 
be related to the temperature exponents as (3 — v{z + 1). 
Here x = 8/E 1 !^ and y — S/T 1 ^ uz are called the scaled 
variables. 

In cases where either the .E-field or the tempera- 
ture has been judged negligible compared to the other, 
data has been examined for one-parameter scaling, i.e. 
whether or not it collapses to a function of only one scaled 
variable. At the core of all arguments for a phase transi- 
tion lies the fact that one in many cases can find a critical 
density n c and exponent, vz or (3, so that data falls onto 
a single curve. This is usually done by trying out differ- 
ent critical densities and exponents. The best fit is then 
determined by eye as the parameter choice that gives the 
best data collapse. Because the question of a phase tran- 
sition in our system still is controversial, we developed an 
objective method to determine the compatibility of the 
data with scaling. 

Our goal is now to give a quantitative measure of the 
consistency of the data with a single curve when we plot 
the resistivities taken at the lowest temperature as a func- 



_ {n. 
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tion of the scaled .E-field, x 

We want to assume as little as possible about the 
form of the scaling function, f(x), but previous inves- 
tigations suggest that Log[f{x)\ is a reasonably smooth 
curve without too much curvature. We therefore investi- 
gate a model where Log[f(x)] is a polynomial with only 
a few terms, f(x) = exp[ao + a\x + aix 1 + a3X 3 + ...]. 
Because theory does not tell us the form of f(x), any 
functional form for / that is general enough is acceptable, 
provided it contains a small number of free parameters 
compared to the number of data points being fitted. If 
we truncate this polynomial at say third order, we ob- 
tain a scaling model for p{n s ,E) that contains the six 
parameters n c , (3, csq, ai, 02 and 03; 



Pmodel(n s ,E) = e «o+aix+a 2 x +a 3X where 

_ (n s - n c )/n c 
x - ^T/p ■ 



(3) 



We can now address the question of the consistency of the 
data with scaling by investigating our ability to fit the 



measured data to this scaling model. We do this with a 
nonlinear regression routine in Mathematica that fits 
all parameters at once and hence obtain values for n c and 
(3 with statistical error estimates. At the same time we 
get an objective measurement of the 'goodness of fit' via 
the x 2 of the fit. 

In order to have an unbiased fit and get x 2 m the cor- 
rect units we need to keep track of the estimated uncer- 
tainty for each experimental point. In this case we need 
to consider not only the error in p itself but also the un- 
certainty in n s and in E. We find that the uncertainty in 
density, a n normally dominates the uncertainty of points 
in the /(x)-plot. To take this into account we propa- 
gate all errors to a single effective error in the resistivity, 

Op.efT 
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dp n 
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dV 



(4) 



We estimate the error in the resistivity a p by propagat- 
ing the errors from the probe voltage ay either via the 
convolution (Q) or the ratio V/I. When calculating the 
total effective error a e g-, we still assume a p and ay to 
be independent. We then weigh the points in the fit with 
a weight 1/er 2 ff . Because the weights depend on the 
fitted polynomial, we have to start with an unweighted 
fit, then propagate the errors and iterate the fit until a 
self-consistent solution is found. 



IV. RESULTS 
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FIG. 3. Scaled data points, p(x) = V/I, from sample one 
and fitted fifth order model polynomial. The fitted polyno- 
mial is only visible where it is not covered by data points. 
Here a low power cut-off of 50pW was used and the fit has 
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= 1. 



P 



(it) 



X 



50 1.0142 ± 0.0006 

100 1.0138 ± 0.0006 

200 1.0134 ± 0.0006 

400 1.0136 ± 0.0009 

500 1.0143 ± 0.0010 



3.83 ±0.03 
3.76 ±0.03 
3.75 ±0.04 
3.81 ±0.05 
3.89 ± 0.06 



1.074 ± 0.009 1.86 

1.083 ±0.009 1.18 
1.094 ± 0.009 0.70 
1.093 ±0.013 0.45 

1.084 ±0.014 0.38 



TABLE I. Results for sample one. Listed are the low-power 
cut-off P c in pW, the critical density n c in units of 10 11 cm -2 , 
the logarithm ao of the critical resistivity in log h/e 2 and the 
average x 2 for the fit. The errors are estimated from the 
statistical deviations of the data from the fitted model only. 



Using either the derivative dV/dl or the ratio V/I 
as the resistivity p, we have used the scaling model de- 
scribed in the previous section to investigate -E-field scal- 
ing in three different samples. Here we used a fifth order 
polynomial ([|) . We also achieved reasonable results with 
lower order polynomials. Because data is taken at a fi- 
nite temperature, (~ 0.2K) we know that the resistivity 
is dominated by the temperature for small enough exci- 
tation -E-fields; hence, we introduce a low power cut-off 
to look for scaling in the Enfield only. 

Using the ratio V/I for the resistivity p we obtain a 
X 2 reasonably close to one. Figure || shows the result 
from sample one, using a low power cut-off P c = 50 pW. 
The critical exponent we f\tvX in this way is comparable 
to previous investigations Jffltj even though the exponent 
and x 2 do depend somewhat on the particular cut-offs 
used. 

We have most data from sample one. It also shows the 
best behavior; x 2 1S about one and the physically inter- 
esting fit parameters, in particular the exponent /3, are 
relatively stable as the low-power cut-off is increased, see 
Table Q. The results in the other samples show stronger 
dependence on the cut-off. Only for quite large low power 
cut-offs (~ 500pW) are the exponents from different sam- 
ples compatible with each other, see Table O for results 



from sample two. We note that at an excitation power of 
1 nW temperature increases in the sample were observed 
via the temperature monitoring RuOj; resistor. Thus, 
data at these high excitations may correspond to the hy^. 
drodynamic regimeEilS discussed in this context earlier .□ 
Sample one shows rather good consistency with the 
scaling model. Here all the fitted parameters are rea- 
sonably independent of the low power cut-off, P c . The 
number of points in a fit varies from 864 for the small- 
est cut-off to 467 for the largest. Actual variations in the 
fitted parameters are consistent with the uncertainty esti- 
mates based on the statistical deviations from the model 
as the low power cut-off is changed. For example, the 
critical density n c is virtually unchanged for all cut-offs. 
This is also reflected in very small statistical uncertain- 
ties for n c in each fit. The exponent (3 (=v (z + 1)) is 
consistent with (3 = 3.8 ±0.05, and the critical resistivity 
e a ° is approximately 2.9/i/e 2 . The mean square devia- 
tion per point, x 2 , however, does vary substantially with 
the cut-off, decreasing as the cut-off is increased. This 
could be understood as an effect of the finite tempera- 
ture the data is taken at; at low powers the temperature 
is not negligible and the model (0) is too crude to give 
a good fit. It is also possible that there are two regimes 
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ao 


x 2 


50 


0.9717 ± 0.0004 


3.14 ±0.02 


2.164 ±0.006 


0.94 


100 


0.9751 ± 0.0004 


3.25 ±0.03 


2.117 ±0.006 


0.63 


200 


0.9780 ± 0.0006 


3.36 ±0.04 


2.078 ± 0.008 


0.46 


400 


0.9818 ± 0.0010 


3.53 ±0.08 


2.028 ±0.013 


0.37 


500 


0.9826 ± 0.0014 


3.69 ±0.13 


2.018 ±0.017 


0.36 



TABLE II. Results for sample two. 



of nonlinearity as mentioned above. The fact that x 2 is 
less than one for the best fits may indicate that we have 
overestimated the statistical uncertainties in the data. 

Sample two shows a larger dependence on the low 
power cut-off, and only at large cut-offs (~ 400pW) is 
the exponent (3 compatible with the results in sample 
one. Clearly the exponent varies much more than the 
statistical error estimates in each of the fits with lower 
cut-offs. Again \ 2 varies in the same way as in sample 
one. At the largest cut-offs, however, both (3 and % 2 are 
close to the values found in sample one. The critical re- 
sistivity is, however, clearly not the same as in sample 
one, e a ° « 7.5/i/e 2 . 

Sample three shows a behavior similar to sample one 
(exponent (3 w 3.7 and e a ° w 3.9) but with much poorer 
statistics. 

In order to gain a further understanding of the data, 
and our fitting procedure, we have also studied the effect 
of introducing a cut-off in the scaled variable x. By keep- 
ing data where |a;| < a; cu t we effectively cut away points 
that are taken far from the critical density. Note that 
the criteria for a specific density to be "near" or "far" 
from the critical density is determined also by the exci- 
tation field; a small field corresponds to a more narrow 
range of densities. Since x depends also on the critical 
density and exponent given by the fit, this cut has to be 
made self-consistently. Ideally one would expect fitted 
parameters to be independent of the actual value of x cu t , 
provided a; cu t is not too large. Figure ^ shows how the 
fitted exponent (3 varies with x cu t in sample one. Each 
point in the figure corresponds to a complete fit using 
both the |a;| < x C nt cut-off and a low power cut-off (P c )- 
The statistical uncertainty in each fit, shown as errorbars 
on the points, are clearly smaller than the actual varia- 
tion in f3. This behavior may be an indication that the 
idea of multiple regimes is appropriate, but further in- 
vestigations are clearly necessary to settle the question. 
The large uncertainties in the points with very small x cu t 
reflect the fact that there are no longer enough data left 
for a good fit. 

We have also investigated the derivative dV/dl for scal- 
ing. In this case we find similar but more varied re- 
sults with generally higher y 2 . The x 2 's determined for 
these fits also depend on the specific Super-Lancos fil- 
ter used (|lj). The derivative scaling also requires the 
introduction of cut-offs for high power data. This was 
done by discarding data outside three times the FWHM 
(Full Width Half Maximum) of the central structure. Re- 
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FIG. 4. Fitted exponent /3 with statistical error estimates 
versus x cut for sample one. Here each point results from one 
fit using a cut-off |x| < a; cu t in the scaled variable x. Points 
with with the same low power cut-off P c are joined as indi- 
cated in the figure. 

suits depend strongly on this FWHM-cut. The results, 
presented above, based on ratios show very little depen- 
dence on such high power cut-offs and we have used all 
data consistent with other cut-offs. 



V. CONCLUSIONS 

We have developed and used a general scaling model 
to fit all parameters in _E-field scaling at once. We find 
some results that are consistent with scaling, \ 2 of about 
one, but also surprisingly large variations in the fitted 
parameters when different physically appropriate cut-offs 
are introduced. 

It would be desirable to generalize and use this method 
for a full two-parameter scaling test, treating the -E-field 
and temperature at equal footing. Any one-parameter 
scaling suffers from the fact that the other (unsealed) 
parameter is not zero. As we have discussed above, any 
Enfield scaling will always suffer from the finite temper- 
ature of the surrounding bath. In the same way, temper- 
ature scaling is sensitive to the nonzero excitation field 
used to measure resistivities, and heating of the 2DES 
will always be a problem for low enough temperatures. 
Both these problems are reduced in a two-parameter scal- 
ing test where the scaling model (0) is supposedly valid 
for any E and T such that a; and y are close enough to 
the transition. 
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